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Abstract 

Closure seems to be something rheologists would prefer to avoid. Here, the story of closure 
is told in such a way that one should enduringly forget any improper undertone of "uncontrolled 
approximation" or "necessary evil" which might arise, for example, in reducing a diffusion equation 
in configuration space to moment equations. In its widest sense, closure is associated with the search 
for self-contained levels of description on which time-evolution equations can be formulated in a 
closed, or autonomous, form. Proper closure requires the identification of the relevant structural 
variables participating in the dominant processes in a system of interest, and closure hence is 
synonymous with focusing on the essence of a problem and consequently with deep understanding. 
The derivation of closed equations may or may not be accompanied by the elimination of fast 
processes in favor of dissipation. As a general requirement, any closed set of evolution equations 
should be thermodynamically admissible. Thermodynamic admissibility comprises much more than 
the second law of thermodynamics, most notably, a clear separation of reversible and irreversible 
effects and a profound geometric structure of the reversible terms as a hallmark of reversibility. 
We discuss some implications of the intimate relationship between nonequilibrium thermodynamics 
and the principles of closure for rheology, and we illustrate the abstract ideas for the rod model of 
liquid crystal polymers, bead-spring models of dilute polymer solutions, and the reptation model 
of melts of entangled linear polymers. 
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I. INTRODUCTION 



The term "closure" usually comes in the combination "closure approximation," and it is 
considered as a necessary evil in deriving an autonomous, or closed, set of evolution equations 
in an attempt to simplify or solve a problem. Employing a closure approximation, usually 
of unclear and uncontrollable quality, then appears as the desperate way out of a dead end. 

The prototype scenario for closure approximations is the passage from probability den- 
sities to moments. The evolution of a probability density, or configurational distribution 



function, of a Markov process is govern ed by some kind of Kolmo g orov 



ential Chapman- Kolmogorov equation Gardinerl (119901 ): lOttinged (119961 )]. say a Fokker 



brward or differ- 



Planck equation (or, a diffusion equation in configuration space). The goal is to de- 
rive autonomous evolution equations for a set of moments. The dilemma is that the 
evolution equations for the moments following from the Fokker-Planck equation involve 
more complicated moments than one intended to consider. Adding these more compli- 
cated moments to the list of variables does not help because their evolution involves 
even more complicated moments. To avoid an intractable infinite hierarchy one needs 
to approximate the more complicated moments in terms of the simpler moments that 
one actually wants to work with. Only such a closure approximation leads to an au- 
tonomous set of evolution equations for a set of moments. In rheology, a number of 
syste ms are treated by clos u re a.pproxi m ations of this type, suc h as liquid-crystal poly- 



mers 



Hinch and Leal 




Larson 



199C): 



(11976h: 



Doil (Il98lh: 



Doi and Edwards! (|l986l) : 



Larson and Ottingerl ( 



Chaubal and Leall ( 



19981 



Feng et al. 



Forest and Wand (|2003|): iKroger et al 



Zubarev and Iskakova 



991) 



Advani and Tucker III 



Bhave et al. 



(Il998f ): 



Ilg et al. 



20081)] , magn etic fluids 



1993)i 



19991 



extensible nonlinear ela s tic ( F ENE) springs 



Wedgewood and Ottingerl (|l9881): 



Martsenyuk et al 



et al 



mi); 



Chauba. 
EdwardJ ( 20021) : 



20001 ): lllg and Kroger! (!2002r)]. or th e nonlin ea r effects of 



Bird et a. 



Wedgewood et al. 



( 199l!): 



1980) 



Ottingei 



Herrchen and Ottingei (|l997 ): 



finitel y 



1987f); 



3i20oJ 



Lielens et al. 



(11999! ): 



Yu et al 



drodynam i c inter ac tions 



ume 



Ottinger! (Il989a): IWedgewoodI ( 



200% 



Ottingei] (Il987 



Prakash and Ottinger! ( 



Du et al 



19891): 



1999) 



112005) 



Ottingei 



Prabhakar and Prakasl 



19896!): JZvlka and Ottingei! (!l989f ): 



Prakash and Ottingei ( 



Prabhakar and Prakash (2002) 



hu 



199?!)] and excluded vol 



Prakash! (!2002!)] in di- 



lute polymer solutions [Bird. Curtiss. Armstrong and Hassager! (!l987! ); 



Bird and Ottingei 
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fll992h]. The 



was 



PeterlinI (11961 



pioneer of closure for a dumbbell model with a nonlinear spring force law 



19661 ). Experience based on a comparison of the solutions of Fokker- 



Planck equations and moment equations suggests that there seems to be enormous room for 
criticism and improvement in the closure game. 

A similar closure problem arises if one starts from a probability density in a high- 
dimensional space and passes to a contracted lower-dimensional distribution by integrating 
out degrees of freedom. A classical challenge is the derivation of evolution equations for the 
single- and two-particle distribution functions from the Liouville equation for the probability 
density of a macroscopically large number of particles. For pairwise interactions, the equa- 
tion for the n-particle distribution function involves also the (n + l)-particle distribution 
function. One thus arrives at the BBGKY hierarchy of equations for contracted distribution 
functions introduced independently by Bogolyubov (1946), Born and Green (1946), Kirk- 
wood (1946), and Yvon (1937). In this context, Boltzmann's celebrated Stofizahlansatz, in 
which the two-particle distribution function at the beginning of a collision is expressed as a 
product of single-particle distribution functions, is the basis for the successful derivation of 
a closed kinetic equation for the single-particle distribution. Also in this context, closure has 
led to doubts and controversy, for example, about its role in the emergence of irreversibility. 

At this point, it should be quite obvious why closure is often perceived as a necessary evil 
and as a questionable but unavoidable mathematical trick with threateningly far-reaching 
consequences. Why then would anybody associate stupendous beauty with closure, as sug- 
gested in the title of this article? As illustrated above, closure has to do with the search 
for simplified autonomous levels of description, and proper simplification is the key to un- 
derstanding by focusing on the essence of a problem. The search for closure, understood 
as the physically motivated search for simple autonomous levels of description for a given 
range of phenomena of interest, is at the heart of developing insightful theories leading to 
fundamental understanding and useful applications in rheology and many other branches 
of science and engineering. I see the beauty of closure in this association with recognizing, 
highlighting and formulating the essentials, as I want to elaborate in this article. I try to 
offer some colorful tesserae which, hopefully, the reader can use to compose an appealing 
mosaic of the role of closure in the process of understanding. 



3 



II. THE BEGINNING: DIFFERENTIAL CONSTITUTIVE EQUATIONS 



The origin of rheology as a scientific discipline lies in the 1920s when Eugene C. Bingham 
coined the term rheology (1920) and the Society of Rheology was founded (1929). Closure, 
in a wide sense, has been a central topic of rheology ever since the derivation of the con- 
vected Maxwell model from molecular considerations by I Green and TobolskyI (119461) and 



the famous formulation of differential rheological equations of state by 



Oldrovdlfll950[ ). Dif- 



ferential constitutive equat ions for the stress tensor have a long and successful history in 



rheology [see, for example. 



Bird. Armstrong and Hassagerl (119871 )]. Such constitutive equa- 



tions successfully describe a variety of nonlinear viscoelastic phenomena in a simple setting, 
and hence they contribute significantly to our understanding in rheology. 

The flow behavior of Newtonian fluids can be described in terms of the five hydrodynamic 
field of mass density, velocity (three components), and temperature (instead of the velocity, 
one can use the momentum density and, instead of the temperature, one can use the internal 
energy or entropy density). This is a natural, universal setting suggested by conservation 
laws. What is the minimum setting for an autonomous description of more complex fluids? 
For rheologists, it seems natural to employ the stress tensor as an important further variable, 
and the search for differential constitutive equations expresses the belief that, with the stress 
tensor alone, one can obtain a closed description of complex fluids. An additional tensor as 
a structural variable provides a classical setting for the autonomous description of complex 
fluids in rheology. 

While we rheologists grew up with differential rheological equations of state, the possi- 
bility of formulating realistic equations of this type is not a priori obvious. The primary 
source of stresses are forces acting over distances, and the stress tenso r can be expressed as 



the av erage of a tensor product of relative position and force vectors Irving and Kirkwood 



( Il950l )]. A differential equation of state implies the possibility of finding a closed description 
for such an average. Describing complex fluids with just an additional tensor always implies 
a closure assumption for an average of some tensor product of two vectors. This becomes 



Bird, Curtiss, Armstrong and Hassaeer 


Peterlin 


1961, 


1966 


) on nonlinear spring 



natural and convenient to work with a conformation tensor, such as the second moment of 
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a dumbbell vector, rather than with the stress tensor directl y. Of course, the stre ss tensor 



Dressier et al. 



fll999h ]. 



must then be expressed in terms of the conformation tensor 

In a wider sense, any choice of an autonomous level of description admitting a closed 
description requires some form of closure assumption. For example, when generalizing from 
dumbbell models of polymer solutions to bead-spring-chain models one looks for closure 
in a whole array of conformation tensors, where Gaussian approximations are particularly 



popular and successful. Models with two coup 
by abstract orientation an d elongation tensors 



fitting of rheological data Edwards et al.l (119961 )]. It is natural to include the heat flux in the 



ed tensor variables have also been motivated 

)] or by the successful 



Wilchinsky et al. 



(12001 



same way as the momentum flux into t he list o : 



tor is hence anot 



ler natural scenario 



(Il996h : 



Jou et al. 



variables. Working with a tensor and a vec 



Miilleil (1196 (t): Beris and Edwardsl (Il994h : 



(Il998h : I Jou and Casas-VazquezI (l200lh : 



Lebon et al. 



Jou et al 



(120081 )] which is use- 
ful not only for comple x fluids u n dergoi ng non-isothermal flow but also to formulate rela- 
tivistic hydrodynamics Ottinged (jl998a Jd)]. The entire framework of extended irreversible 
thermodynamics has been buil t on the idea to e mploy the morn er itum and heat fluxes as 
additi o nal structura l varia bles |Jou et al. ( 1996 ): Jou et al.i (il998i) : Ijou and Casas- Vazquez 



feoOll); 



Lebon et al 



eral form in 



Ottingei 



( 2008h 



The level of a tensor and a scalar has been analyzed in its gen- 



(I2OO2I ) because it has found a nur nber of modern app 



ample , in the Doi-Ohta model of emulsions and blends Doi and Ohtal (Il99ll ): 



ic ations, for ex- 



Wagner et al 



(1l999l )]. where the scalar describes the amount of interface per unit volume and the tensor 
variable separately accounts for the orientation of the in t erface, or i n the pompon model for 



melts of branched polymers [McLeish and LarsonI (jl998|); lOttingerl (l200ll )]. where the scalar 
represents the stretch of a tube confining the molecular backbone and the tensor is used to 
characterize the tube orientation, or in transient network models of polymer melts, where 
the scalar describes the number density of network segments and the tensor characterizes the 
segment deforr nation. An even s impler example is prov i ded b y inhomogeneous dilute poly- 
mer solutions Ottingerl (1l992l ): iBeris and Mavrantzad (1l994l )]. where the scalar describes 
the polymer concentration and the tensor represents the polymer stretch and orientation. 

In all the above-mentioned theories with a tensor variable, underlying closure assumptions 
must exist, at least, if one models nonlinear effects. The success of differential equations of 
state suggests not to consider closure as a necessary evil but as an important step toward 
grasping the qualitative essence of rheological phenomena. If one thinks about the Doi-Ohta 
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or pompon models one appreciates the understanding of the rheological behavior of complex 
systems based on the evolution of structural variables, and one does not at all worry about 
closure approximations. As a matter of fact, closure creates comprehension. 



III. THERMODYNAMIC STRUCTURE 

To understand the essence of rheological and other nonequilibrium properties, one should 
formulate equations on as coarse as possible levels of description. Passing to a coarser level of 
description should not be regarded as a formal mathematical approximation but rather as an 
insightful physical identification of the relevant variables, that is, as an intellectual achieve- 
ment. We are then faced with the following question: What are the fundamental physical 
principles that should be respected in formulating an evolution equation on any autonomous 
level of description? This is, of course, the realm of nonequilibrium thermodynamics. In the 
ideal case, the fundamental thermodynamic structure enforces not only physically admissi- 
ble but also mathematically well-behaved equations for which the existence and uniqueness 
of solutions can be demonstrated. The thermodynamic structure should also be preserved 
in numerical integration schemes. 

After respecting all conservation laws, the most prominent thermodynamic principle is 
certainly the second law excluding the possibility of a negative entropy production rate. 
Also the clear separation of reversible and irreversible contributions to evolution equations 
is of fundamental importance. The reversible contribution should possess a rich structure 
reflecting the idea of "mechanistic control," and it should not touch the entropy. The 
reversible contribution is generally assumed to be of the Hamiltonian form and hence requires 
an underlying geometric structure (given in terms of a Poisson bracket or Poisson operator). 
The remaining irreversible contribution is driven by the nonequilibrium entropy by means of 
a dissipative bracket. In the GENERIC ("general equation for the non equilibrium reversible 



i rreve rs ible coupling" ) framework o 



n onequilib r ium t hermo dynamics 



(I1997I ): lOttinger and Grmelal (119971 ): lOttingerl (120051 )]. these ideas are condensed into the 



Grmela and Ottinger 



evolution equation 

dx L ^'^ M (V\ 
dt Sx Sx'' 

where x represents the set of independent variables required for a complete description of a 
given nonequilibrium system, E and S are the total energy and entropy expressed in terms 
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of the system variables x, and L and M are certain linear operators, or matrices, which 
can also depend on x. Equation ([1]) is supplemented by the complementary degeneracy 
requirements 

and 

M.f.O. (3, 

The requirement that the entropy gradient 6S/6x is in the null-space of L in Eq. ([2]) expresses 
the reversible nature of the L-contribution to the dynamics: the functional form of the en- 
tropy is such that it cannot be affected by the operator generating the reversible dynamics. 
The requirement that the energy gradient 6E/6x is in the null-space of M in Eq. ([3]) ex- 
presses the conservation of the total energy in a closed system by the M-contribution to the 
dynamics. The two contributions to the time-evolution of x generated by the energy E and 
the entropy S in Eq. ([T]) are called the reversible and irreversible contributions, respectively. 

Further general properties of the matrices L and M are discussed most conveniently in 
terms of the Poisson and dissipative brackets 

l-4.Bl = |^-M.^, (5) 

OX ox 

where A, B are sufficiently regular real-valued functions on the space of independent vari- 
ables. In terms of these brackets, Eq. ([1]) and the chain rule lead to the following time- 
evolution equation of an arbitrary function A in terms of the two separate generators E and 

HA 

— = {A,E} + [A,S]. (6) 
The further conditions for L can now be stated as the antisymmetry property 

{A,B} = -{B,A}, (7) 

and the Jacobi identity 

{A, {B, C}} + {B, {C, A}} + {C, {A, B}} = 0, (8) 
whereas the product or Leibniz rule for Poisson brackets, 

{AB,C} = A{B,C} + B{A,C}, (9) 



follows immediately from the definition in Eq. (jl]). In these equations, C is another arbitrary 
sufficiently regular real-valued function on the state space. These properties are well-known 
from the Poisson brackets of classical mechanics, and they capture the essence of reversible 
dynamics. The Jacobi identity ([H]), which is a highly restrictive condition for formulating 
proper reversible dynamics, expresses the invariance of Poisson brackets in the course of 
time {time- structure invariance). 

Further properties of M can be formulated in terms of the symmetry condition 

[A,B] = [B,Al (10) 

and the non-negativeness condition 

[A,A]>0. (11) 

This non-negativeness condition, together with the degeneracy requirement ([2]), guarantees 
that the entropy is a nondecreasing function of time. 

The properties fllOl) and (ITT]) imply the symmetry and the positive-semidefiniteness of M 
[for a more sophisticated di scussion o f the O nsaqer-Casimir symmetry properties of M, see 



Sections 3.2.1 and 7.2.4 of 



Ottingerl (l2005l )]. From a physical point of view, M may be 



regarded friction matrix. 

The thermodynamic structure summarized here can also be regarded as a geometric struc- 



ture. In the mathematical 

metriplectic 



Morrison 



iterature, this geometric structure is sometimes referred to as 



(119861 )]. In particular, Poisson operators are in ti mately related to (du 



Marsden and Ratiu 



als of) Lie algebras [see, for example. Appendix B of lOttingerl ([20051) or 
(119991 )]. For example, convection effects are related to the Lie group of space transfor- 
mations and its representations; as a result, we obtain convected rather than partial time 



derivatives, as required by the 



Bird. Armstrong and Hassagerl (119871 ): iLodgd (Il974j )] 



a mous prmcip 



e of m aterial objectivity or frame indifference 



IV. DOI AND BINGHAM CLOSURES 

We now have the thermodynamic tools to address the topic of closure thoroughly. One 
of the most famous closures in rheology is the quadratic ansatz 

— 7 : {uuuu) = —7 : {uu){uu) (13) 
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introduced by Doi in the t heory of hqu i d crystal p olymers Doil (jl98ll )]. The Hess-Doi theory 



of liquid crystal polymers 



HessI fll976h 



Doi 



(Il98ll )] is based on a configurational distribution 
function f{u) for an ensemble of rigid rods, where u is the orientation vector of a rod. The 
averages in Eq. f[T^ are performed with the configurational distribution function f{u). If 
7 is the sum of the velocity gradient tensor and its transpose, the left-hand side of Eq. (IT5]1 
constitutes a convective contribution to the time evolution of the second moment tensor 
(uu). As one wishes to obtain a closed evolution equation for (uu), the fourth moment has 
to be expressed in terms of second moments, and Eq. (fT3|) offers the simplest possibility. An 



appealing alternative is to assume 



{uuuu) = J uuuu f{uu) (w) d' 



(14) 



where //„.„.) (u) is a given cl ass of distribution functions parametrized by the second moments 
Chaubal and Leal ( 1998 ): Feng et al] (1998)]. The most popular choice is the exponential 
of a q uadratic form of u, the mathematical-statistical properties of which have been studied 



by C. 



Bingham! (119741 ) (who should not be confused with the pioneering rheologist Eugene 



C. Bingham). 

The idea of introducing parametric distributions to obtain closure is widely used in 
nonequilibrium statistical thermodynamics. Generalized canonical distribution functions 



paramet rized by 



straints 



Ilg et al 



agrari g e multipli e rs are obtained by maximizing the entropy under con- 



(12002); 



Ilg et al 



(120031 )]. For example, th e Bingham dist ribution arises 
by maximizing the entropy for a fixed second-moment tensor |Ilg et al.i (|2003l )]. The deriva- 
tive of the entropy with respect to the second-moment tensor is then given by k-^A, where 
A is the Lagrange multiplier associated with (uu) (the Bingham distribution is propor- 
tional to exp{— A : uu}). Efficient integration schemes and accuracy control, with the 
possibility of changing the level of description upon a loss of acc uracy, ha v e been proposed 



i n the manifold of generalized canonical distribution functions Ilg et al.l (120021 ): 



Ilg et al. 



( 120031 )]. which is also known as the quasi-equilibrium manifold. The idea of parametric 



distributions has been further develop e d into the powerfu 



fold method 



Gorban and KarlinI (119921 ): 



Gorban and Karlin 



tool of the invariant mani- 



(!l994h 



Gorban et al. 



mm; 



Gorban and KarlinI ((20051)], where geometric ideas and thermodynamic projectors are found 
to offer a more elegant and general approach to closure problems than explicit parametriza- 
tions, and an iterative Newton method conveniently provides successive improvements of 
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the equations for the moments. However, the invariant manifold method cannot produce 
any dissipation; in particular, it cannot produce irreversible equations from reversible ones. 
To allow for this po ss ibility, the method has be en enhanced by Ehrenfest coarse graining 



Gorban et al 



fl200lh : iGorban and KarlinI (120051 )]. 



The closure approximations fll3p and fll4l) . among many others, ha ve been studied 



i ntens el y in the 



1976) 



Doi 



19811 



it erature on the rod model of liquid crystal polymers 



Doi and Edwards ( 



1986) 



Advani and Tucker 



linch and Leal 



199G) 



Larson ( 



199G) 



Larson and Ottinge 



f|l998h: 



Feng et al. 



199l|); 



Bhave et al 



(11998h : lllg et al. 



1993); 



(jl999|) 



Chaubal et al. 



1995); 



Chaubal and Leal 



Edwardsl (12002); Forest and WangI (120031); 



Kroger et al.l (j2008l )]. Some of these closures are extremely successful in particular flow 
situations but fail terribly in others. For example, in simple shear flow, Doi's quadratic 
closure f|T3|) . when applied to the fourth moments occurring both in t he reversible and i n 



the irreversible terms, exhibits only time-independent stable solutions 



Chaubal et al. 



Bhave et al. 



(119931 ): 



(119951 )] and thus misses the well-known periodic solutions known as "tum- 
bling" (preferred axis of alignment rotates in the plane of shear) and "wagging" (preferred 
axis oscillates back and forth in the plane of shear), whereas the Bingham closure admits 
such time-dependent stable solutions. Can one identify the ultimate winners and losers by 
such observations so easily? 

At this point, we should remember our goals and ambitions. We are not really inter- 
ested in the Olympic Games of mathematical closure approximations, in which gold, silver, 
and bronze medals are given away in a myriad of specialized flow disciplines. We rather 
strive after fundamental understanding by identifying and verifying autonomous levels of 
description bringing out the essence of a problem, without necessarily reproducing all the 
details. In the present case, the autonomous level of description for liquid crystal polymers 
is proposed to be given by the structural variable (uu), and we must respect the structure 
of thermodynamically admissible equations on this level. Moreover, judging the success of 
closure approximations by a comparison with the exact results for all kinds of flow situations 
is not particularly useful; we should clearly prefer to find a priori criteria for overall success 
that do not require any knowledge of exact solutions. 

As the fourth moment in Eq. (fT3l) is associated with convection, we deal with a reversible 
term generated by the energy gradient with the help of a Poisson operator. As the kinetic 
energy, which generates convection, is unaffected by the structural variable, we need to 
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focus on the Poisson operator and its properties as the hallmarks of reversible motion under 
mechanistic control. The most restrictive criterion is the Jacobi identity ([8]) expressing the 
time-structure invariance of reversible dynamics. This crit erion has been analyzed in great 



detail in a similar context by 



Edwards and Ottingerl (119971 ). and we here merely summarize 



and discuss the r nost important results. 



According to Edwards and Ottingen (Il997l ). reversible dynamics obtained under the 
Doi closure f|T3|l can be generated by a valid Poisson operator, whereas this is impos- 
sible for the Bingharn closure (IT^ and many other common closures. The analysis of 



Edwards and Ottinged ( 119971 ) proceeds according to symmetry. Note that the fourth mo- 



ment is symmetric in all four tensor indices. This full symmetry, which is respected by the 
Bingham closure, is incompatible with time-structure invariance. The Doi closure possesses 
symmetry in pairs of indices, and under exchanging the pairs. For this lowered level of 
symmetry, the Doi closure (fT3l) is the only possible closure compatible with time-structure 
invariance. For even lower leve ls of symmetry, further admissible closures can be constructed 
Edwards and Ottingerl ( 119971 )]. 

Of course, you could say "Why should I, as a rheologist, make my life more complicated 
by worrying about some mysterious Jacobi identity?" Because Poisson structures are at the 
heart of reversible dynamics, and because violating laws out of ignorance does not protect 
you from punishment! Of course, one could still try to argue why a law is not applicable 
in a particular situation, but the only reason for losing mechanistic control in reversible 
dynamics that I am currently aware of is the presence of nonholonomic constraints. 

If one has to choose or compromise between symmetry and time- structure invariance, why 
should symmetry be less important? The answer should be clear by now: Because our goal is 
to establish a healthy autonomous level of description based on second moments, consistent 
with all the laws of nonequilibrium thermodynamics, and not necessarily to achieve a faithful 
mathematical approximation of a fourth moment in terms of second moments. 

If only the quadratic closure (JT3l) is thermodynamically admissible, does thermodynamics 

i ng" a n d "wagging" in 



thus force us into an inferior closure that cannot even p redict "tumb^ 



shear fl ow? A much more b alanc ed vie w was offere d by 



Bhaveetal 



( 1l995l ). IChaubal and Leall (119981 ) . and iFeng et al.l (119981 ) . based on a "solution map" or 



(Il993f ). 



Chaubal et al 



"bifurcation set" in the parameter space of two-dimensional flows and nematic strength. 
Shear flow appears as a singular special case, and "tumbhng" and "wagging" do actually 
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occur within the quadratic closure for flows that are only very slightly more rotational than 
simple shear flow. The overall solution map is deformed only slightly for the quadratic 
closure, but with rather serious consequences for the singular special case of shear flow. 
Moreover, none of the known closures exhibiting "tumbling" and "wagging" can predict 
the proper t ransit ion from "wagging" to flow aligned steady solutions at high shear rates 



Feng et all (119981 )]. 



A more rotational behavior with "tumbling" and "wagging" can also be achieved by going 
from the upper convected codeformational derivatives appearing naturall y in the second mo- 



Beris and Edwards 



ment equations to mixed or Schowalter derivative s [see pp. 556 an d 568 of 
(|l99J)]. It has been elaborated in Section 4.2.3 of lottingerl ko03^ that, within the GENERIC 
framework, such a modiflcation can be implemented through an additional irreversible con- 
tribution, as suggested by the occurrence of a slip coefficient. Schowalter derivatives are 
represented by an antisymmetric contribution to the friction matrix so that, according to 
Eq. ( fT2l) . they do not lead to entropy production. 

We have focused entirely on the fourth moment (|T3il that appears as a convective con- 
tribution to the time evolution of the second moment tensor (uu). We have not paid any 
attention to the fact that, for the Maier-Saupe mean-fleld nematic potential, there occurs 
another fourth moment in the irreversible contribution to that evolution equation so that a 
further term requires closure. As a matter of fact, a different closure could be used there 
because irreversible contributions are much less restricted than reversible ones. For example, 
one could use the Doi closure in the reversible contribution to fulflll time-structure invari- 
ance, and the more symmetric Bingham cl osure in the i rrever sible contribution in order 
to combine the advantages of both closures [Sgalari et al.l (120021 )]. where the less restricted 
closure in the irreversible term seems to be the bigger source of problems. Actually, the fact 
that the GENERIC framework of thermodynamics expresses time evolution in terms of the 
generators E, S and the matrices L, M strongly suggests to consider each of these building 
blocks separately, that is, with the natural possibility of separate closures in the reversible 
and irreversible terms. In the subsequent section, we look at the implications of thermo- 
dynamics for closure in the irreversible contribution to the time evolution of moments. We 
do this in the context of the Gaussian approximation, which has found many successful 
applications in polymer kinetic theory. The rod model of liquid crystal polymers could be 
investigated in a similar way, but the handling of Gaussian distributions is more familiar. 
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V. GAUSSIAN CLOSURE 



In the Gaussian closure procedure, or Gaussian approximation, the Jacobi identity is 
not an issue. Because there are no constraints and Gaussian objects are deformable, one 
can reahze codeformational behavior in the forr n of upper conve cted derivatives, which are 



known to be consistent with a Poisson bracket [Ottingerl (120051 )]. Nevertheless, something 
remains to be checked for the irreversible term. As Gaussian approximations are usually 
implemented on the level of time-evolution equations, one needs to check whether they are 
consistent with a Gaussian entropy generating irreversible dynamics via the friction matrix 
obtained for Gaussian distributions according to Eq. ([T]). 

To discuss the proper formulation of the irreversible term in the Gaussian approximation, 
we start from a Fokker-Planck equation of the general form 



df d f d 



dt dxj \ 2 J 2 dxj dxk ' 

where / = f{x) is a probability density in some i^-dimensional space with coordinates 
Xj, Aj is a reversible drift vector, the positive-semidefinite symmetric diffusion matrix Dj^ 
describes the irreversible effects, 

d In f ^1 

is the effective force implied by the equilibrium probability density f^^, and Einstein's sum- 
mation convention is assumed (one needs to sum from 1 to over all indices occurring 
twice) . We look for (because we believe in the adequacy of) a description on the level of a 
matrix c of second moments with entries, 

Cij {^XiXj^ j XiX j f d X, (-^'^) 

where the coordinates Xj are assumed to be introduced such that = 6ij. We further 
assume that the first moments vanish for symmetry reasons. In the Gaussian approximation, 
these assumptions imply Fj = xj. 

The entropy on the level of second moments can be obtained by evaluating the Boltzmann- 
type conformational entropy —k^ J f \n{f / f'^'^)d^x, where is Bolt zmann's constan t, for 



Gaussian distributions. The result is [see, for example. Exercise 66 of lOttingerl (120051 )]: 



S = ^kB {K - Cjj + In det c) , (18) 
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with the derivative 



dS 

dcij 



(19) 



Equation ( |T8l) for the entropy is an essential feature of any Gaussian approximation. 

For the exphcit transformation rule (fT7j) from probability densities / to second moments 
(xiXj), one can transform the friction matrix M occurring in the irreversible contribution to 
the Fokker-Planck equation f|T5l) according to Eq. (6.180) of lOttingerl (120051 ). The result is: 



2kBMij^ki = (xiXkDji) + (xiXiDjk) + {xjXkDa) + {xjXiDik). 



(20) 



All four contributions on the right-hand side of this equation are equivalent if Mij^^i is con- 
tracted with symmetric tensors; otherwise, the four contributions imply that the contraction 
should be done with symmetrized tensors only. The Gaussian approximation can now be 
introduced into Eq. fl20p b y using Wick's theorem to reduce the order of the moments (see, 
for example, Eq. (2.61) of lOttingeii (119961 )). We successively find 



and 



{XiXkDji) = Cik{Dji) + Cmk ( X, 



{XiXj^D ji^ Ciki^D jij -\- CmiCnk 



dXra I 



(21) 



(22) 

represents the 



^ dxijidx^i I 

The first term on the right-hand side of each of the Eqs. (I?I1) and 
effect of a self-consistently averaged diffusion matrix, whereas the second term accounts for 
fluctuation effects. The description of fluctuation effects by second-order derivatives looks 
particularly natural. Equation (|2T]) can be rewritten in the alternative form 

dDjk \ 



Cj^l{xiXkDji) = {Dij) + l^Xi 



dxk 



(23) 



As we have found the natural entropy and friction matrix for any Gaussian closure, 
we can now compare to the evolution of the second moments obtained directly from the 
Fokker-Planck equation ([T^ . 



^ = {^i^i) + {xi^j) - ^{DikFkXj) - ^{xiFkDkj) + {Dij) + ^ Ui^^ 



1/ dD,k 

- Xi 



dxk 



(24) 



With Wick's theorem in the form of Eq. (l23l) and the Gaussian property Fj = Xj, Eq. ([24 
can be rewritten as 



= {AiXj) + {xiAj) + -{{xiXkDji) + {xjXkDu)) (c^/ - Ski) ■ 



(25) 
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The irreversible contribution to this evolution equation is exactly what one recovers from 
the irreversible contribution to GENERIC by combining Eqs. f|T9l) and fl20l) . The Gaus- 
sian approximation is thus nicely consistent with an irreversible contribution to dynamics 
generated by the Gaussian entropy on the level of second moments. Contrary to our dis- 
appointing experience with the reversible term, parametric density estimation works very 
nicely for Gaussian approximations to the irreversible term. 

In the Gaussian approximation, a configurational distribution function / is assumed to be 
Gaussian at any time t. This does not necessarily imply a Gaussian stochastic process, for 
which all joint distributions at different times must also be Gaussian. The construction of a 
full Gaussian process governed by a linear s tochastic differe ntial equ ation has been described 



and discussed critically in Section 4.2.4 of 



Ottingerl (119961 ) and by 



Hiitter et al 



fl2003h . A 



well-defined and physically consistent stochastic process on the level of second moments can 
)e introduced by ad ding noise to Eq. (!25|) according to the fluctuation-dissipation theorem 
Hiitter et all (2003)]. 

Because the Gaussian approximation has been used very successfully in the kinetic 
theory of dilute po l ymer s ol utions, where it h a s bee n applied t o hyd ro dynamic in 



teractions 



198911: 



Ottingei ( 



1989i) 



Prakash and Ottinger 



Zvlka and Ottingeil fll989f ): lOttingei] f 



Prabhakar and PrakashI (120021 ): 



19971)]. excluded volume 



1989fl); 



Wedgewood 



Prakash and Ottingei ( 



PrakashI (120021 )]. and internal viscosity 



Schiebei 



19991); 



we here specialize our gene ral results to t he ca se of bead-spring chains, which has pre- 



viously been considered by 



Hiitter et al. 



(120031 ). If the chains consist of beads, we 
have K = 3{N — 1) configurational degrees of freedom and the matrix c is an array of 
(A^ — 1) X (A^ — 1) tensors representing the variances and covariances of the A^ — 1 
connector vectors, where the normalization condition at equilibrium now reads c^^ = 6jkl- 
According to Eq. (fTSll . the conformational entropy per polymer molecule is given by 



2 



N-l 



tr(l — Cjj) + In det c 



(26) 



where c is the large matrix consisting of (A^ — 1) x (A^ — 1) blocks Cjk, each of which is 
represented by a 3 x 3 matrix. As the GENERIC framework involves only derivatives of 
entropy, we provide the result ([TI?]) in the form 

dsp 

dcjk 2' 



k - ^jfci 



(27) 
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where the tensors cL are represented by the blocks of the inverse of the large matrix c, 



N-l 



(28) 



1=1 



The degeneracy condition ([2]) of nonequilibrium thermodynamics implies that the entropic 
spring contribution to the pressure tensor is given by 



N-l 

U = npT{ Spl + J2 

j,k=i 



' dsp 



+ [Cjk 



■jk 



(29) 



where T is the absolute temperature and rip the number density of polymers. By inserting 
the derivatives (127|) into Eq. (!29|) and using the symmetry of the large matrix c (which 
implies (cj^)^ = c^j and an analogous identity for the inverse), we obtain 



N-l 



n = HpSpT 1 + UphsT '^n)- 



(30) 



Note that this simple form of the pressure tensor consisting of Hookean spring contributions 
from each of the connectors is a direct consequence of the functional form of the entropy in 
Eq. ([26]). 



VI. VARIABLES ARE EVERYTHING 



For the development and discussion of the Doi, Bingham, and Gaussian closures in the 
preceding sections, we looked in great detail at complicated moments and their expression 
in terms of second moment tensors. Over all these details we should not forget that our sole 
goal was to establish an autonomous level of description based on second moment tensors. 
The important questions are: Can this be done? How can this be done? An important 
outcome of the general analysis of the second-moment level of description is: There is not 
much choice in the reversible dynamics. 

The key problem of nonequilibrium thermodynamics or, if you like, of coarse graining, 
or of closure, or of understanding, is the choice of good variables for a problem of interest. 
For this crucial task of choosing variables, thermodynamics leaves you alone with your 
insight, intuition, imagination, and ingenuity (i'*). To gain insight, of course, experimental 
results are of particular importance. Once you have chosen your variables and expressed 
them in terms of the variables of a more detailed well-established level of description, often 
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the atomistic level, statistical thermodynamics p rovides systemat ic 



recipes 



'o r calc u lating the 



therm o dynamic buildiri q blocks E, S, L, and M |Ottingerl (Il998d ): lOttingerl (|2000| ): lOttingei 
( I2OO5I ): lOttingeii (I2OO7I ): lllg et al.l (120091 )] . These building blocks imply the autonomous time 
evolution ([T]) and hence all closure properties. The generalized microcanonical, canonical, 
and mixed ensembles of nonequilibrium statistical mechanics can be considered as natural 
candidates for parametric density estimation. Any concrete realization of the GENERIC 
structure provides a solution to the closure problem. 

The importance of the choice of variables can be illustrated nicely for the reptation 
model of melts of entangled linear polymers. One usually does not speak about closure in 
this context, however, one clearly uses mean-field type and further simplifying assumptions. 
Starting from the picture of a large number of massively entangled random walk chains, one 
first assumes that one can look at the motion of a single probe chain constrained by a tube 
or slip links produced by other chains, or anisotropic friction accounting for the hindrance of 
sideway motions by other chains. In a second step, the single-chain picture is further reduced 
to that of a single segment with orientation it at a position s within the chain, where the 
label s varies from to 1 in going from one chain end to the other. In order to achieve this 
further simplificatio n, one needs to make an as s umption like "independeii t alignment of large 



Doi and Edwards ( 



straigh t segments" 



( ll978d )] or "smooth curvature" 



1978fl), 



Doi and Edwards ( 



Curtiss and BirdI ( 



imJ) 



19786). 



Doi and Edwards 



Curtiss and BirdI ( 



198161 )]. The 



entanglement length scale is thus introduced as the length of the independently aligned 
segments or as the persistence length of smoothly curved chains. The truly heroic idea 
of the reptation model is to postulate that the single- segment configurational distribution 
function f{u, s) leads to an autonomous level of description for the complicated system of 
entangled chains. This postulate implies a natural but highly nontrivial closure. 

The formulation of the thermodynamic building blocks, and hence of the evolution equa- 
tions, is a much simpler step than the choice of variables and can actua l ly be done with very 
little of i^, as has been shown by the author in Section 8.4.6 of lOttingerl (l2005l ). Once f{u, s) 
has been identified as a good structural variable, there is so much guidance from nonequi- 
librium thermodynamics that, for example, one is automatically reminded to consider con- 
straint release associated with the reptation of constraining chains and the possibility of 
anisotropic tube cross sections. 
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VII. REDUCTION VERSUS COARSE GRAINING 



For diffusion equations with Gaussian solutions, which are associated with hnear stochas- 
tic differential equations, the second-moment equations reproduce all the features of the ex- 
act solutions. Closure is usually considered as an approach to achieve approximate solutions 
to nonlinear problems, with the ambition to be as accurate as possible. This is the idea 
of a solution or reduction technique. In the context of the invariant manifold method, we 
have already encountered the hallmark of reduction techniques: they do not produce any 
additional dissipation and, in particular, reduction techniques cannot lead from reversible 
to irreversible equations. Equation (I2U]) for the Gaussian approximation also exhibits the 
hallmark of reduction: the friction matrix on the level of moments is directly proportional 
to the diffusion tensor Djk on the level of the configurational distribution function. Also dy- 
namic renormaliza tion, as carried out in the context of hydro dynamic interactions in dilute 



polymer solutions 



Ottingerl (120091 )] . turns out to be a reduction technique. 



To illustrate the closure problem in the introduction, we had also mentioned the BBGKY 
hierarchy and Boltzmann's kinetic equation, where the former is reversible and the latter 
is irreversible. Boltzmann's derivation of his irreversible kinetic equation for rarefied gases 
from the reversible equations of classical mec hanics is an enormous achievement that has 



created a lot of controversy and deep insights 



Lanford III! (Il975l )]. The derivation of Boltz- 



mann's kinetic equation cannot be achieved by a reduction technique. When new dissipative 
processes arise in the passage from a more detailed to a less detailed level of description, we 
speak of a coarse- graining technique. 

In the case of Boltzmann's kinetic equation, the clear separation of two time scales, 
namely the duration of collisions and the time between collisions, is at the origin of irre- 
versibility. A detailed analysis of two-particle collisions is required to obtain the Boltzmann 
equation which itself cannot resolve any processes on the short scale of the duration of a 
collision. In general, in coarse graining one treats fast processes as fluctuations and fluc- 
tuations are associated with dissipation. Projection operators provide a powerful tool to 
separate f ast an d slow process es , thus provid in g the statis t ical m e chanics of coar se graining 



ZwanzieJ fll96lh : iMoril (ll965of ): iMoril (11965JD; iRobertsonI (119661 ): iGrabertI (Il982h : lOttingei 
J2OO5I )]. 

Note that also the diffusion equations for configurational distribution functions used 
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in polymer kinetic theory themselves are the result of coarse graining. In particular, 
spring forces between refer ence points in linear polymer molecules have been discussed by 
Underhill and Doyld (120041 ) by means of statistical mechanics in constant extension and con- 
stant force ensembles. If, however, one is interested in spring forces between beads, each 
located at the ce nter of mass of th e smaller units c oarse grained into the bead, the force 

Jl99ll ). bttingej feooof )]. 



laws are different LasoetaL 



Both reduction and coarse graining lead to autonomous evolution equations as a result 
of passing from a more to a less detailed level of description. Both techniques must preserve 
the GENERIC structure of nonequilibrium thermodynamics, both techniques are associated 
with closure. Whereas the approximate character of closure may be viewed as a disadvantage 
in reduction, it is the actual goal of coarse graining: fast processes are eliminated in favor of 
dissipative processes. This is understanding by focusing on the relevant slow features. One 
should always try to obtain more transparent equations by coarse graining and only when 
further coarse graining is impossible, reduction and explicit solution techniques should be 
used. With the tremendous increase of computer power, there is the temptation of producing 
numerical solutions without exploiting the full potential of coarse graining. Giving in to this 
temptation, we miss the opportunity of deeper understanding and, in the long run, our 
problem-solving skills will deteriorate. 

VIII. RHEOLOGY WITHOUT CLOSURE 



We have recognized the closure problem, in a wide sense, as the search for autonomous 
levels of description on which a closed formulation of the evolution of systems is possible. 
An autonomous system remembers its history, in particular, its flow history, only through 
the current state of the structural variables of the level of description. There is no explicit 
memory and the GENERIC evolution equations ([T]) are pure first-order differential equations. 

In rheology, however, al so integral rheological constitut e equa tions and memory integral 
expansions are well-known Bird. Armstrong and Hassagerl (119871 )]. In view of the enormous 
possibilities of formulating memory functionals, it is very difficult both to select an appro- 
priate memory functional for a given proble m and to formulate the p roper thermodynamic 



principles for memory functionals in general Coleman and Nolll (119631 )]. In particular, there 



is no direct scheme for calculating memory functionals from atomistic models by means of 
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statistical mechanics. Whereas some estabhshed classes of integral models can certainly de- 
scribe rheological measurements for a variety of systems qualitatively well, and there usually 
are sufficiently many fitting parameters to get quantitative agreement, there are no system- 
atic ways of translating molecular understanding of a system into successful memory integral 
models. The identification of key processes and relevant variables, supported by statistical 
mechanics in evaluating thermodynamic building blocks, is a much more promising pathway 
to understanding. 

Of course, there are connections between differential and integral rheological equations of 
state or, more generally, between autonomous levels of description and memory functionals. 
By solving or integrating the evolution equations for the structural variables, in principle, 
one can determine memory functionals. However, explicit solutions can only be obtained 
in exceptional (linear) cases. For example, for the reptation model, a deformation measure 
and a memory function can be identified such that its solution can be written as a simple 
memory integral; the reptation model of Doi and Edwa rds has actually been recognized as 
an integral model of the K-BKZ type [see, for example, iHassagerl ( 1l98ll )]. 

Another famous class of evolution equations with memory are mod e coup li ng th eories. 



As the memory effects are included by following well-defined recipes 



Gotzel fll984h ]. one 



might like to speak of mode coupling as "quasi-closure." In particular, mod e coup 



ory has been foun d to b e useful in the discussion of the glass transition [Gotzd (Il984l ): 



i ng the - 



Kob and Andersen 



(ll99J); 



Kob and Andersen 



(I1995I )]. The glass transition problem is con- 



sidered to be a very difficult one because the ingenious structural variable for describing it 
by autonomous evolution equations has not yet been discovered. The potential energy land- 
scape a nd the i n heren t structure e nergy a re clearly impo r tant c oncepts Stillinger and Weber 



([19821); lAngelll ([19951); IStillinged (119951 ): 



La Nave et al. 



20031 )]. and the inherent structure 
pair correlation function has been proposed as a pr omising c a ndida te for a coarse grained 
description of the approach to the glass transition [Ottingeii (120061)]. This idea has act u- 
ally revealed a static signature of the onset of the glass transition [Del Gado et al.l ([20081)]. 
Structural variables associat e d wit h "shear transformation zones" may be an interesting 



alternative 



Falk and Langeil (119981 )]. 
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IX. CLOSURE 



We have offered some reffections on the topic of closure which strongly suggest that clo- 
sure should not be considered as a helpful mathematical approximation but as the corner- 
stone of establishing autonomous levels of description. Choosing good structural variables 
is the physical prerequisite to achieve closure and to perform successful coarse graining, 
that is, to elaborate the essential features of problems. Closed simplified descriptions lead 
to understanding. We have illustrated this situation in the context of the highly restricted 
fourth-moment closure in the convection term of the rod model of liquid crystal polymers, 
of the admissibility of Gaussian closures in the irreversible term for bead-spring models of 
dilute polymer solutions, and of the far-reaching ingenuity of the variables in the reptation 
model of melts of entangled linear polymers. Liquid crystals and Gaussian closures provide 
examples of reduction, that is, approximate solution procedures. Even more revealing is clo- 
sure obtained from a coarse graining procedure, which replaces fast processes by dissipative 
ones, such as the reptation process. 

With mild exaggeration (rooted in enthusiasm rather than one-sidedness), our reflections 
on closure can be expanded into a full world-map of rheology. The milestones of theoretical 
rheology are ground-breaking new models, and these come with and actually via the insight- 
ful and revealing discovery of good variables. Then, closed equations can be formulated on 
the corresponding autonomous level of description. Nonequilibrium thermodynamics pro- 
vides the structure of admissible evolution equations, nonequilibrium statistical mechanics 
offers the recipes for finding concrete realizations of the thermodynamic structure for given 
variables. We thus obtain all rheological information for the system defined through the 
ingeniously chosen variables. 

Internal energy and entropy are key concepts in nonequilibrium thermodynamics. This 
clearly indicates the way to go for rheology. Focusing on the momentum balance and mo- 
mentum flux has been a good starting point, but a detailed investigation of the associ- 
ated entropy flux must be the next step. This also requires consideration of compress- 
ibility effects. How do we need to generalize the concepts of thermal conductivity and 
heat capacity in complex fluids undergoing flow? Pioneering exper iments on anisotropi c 



heat flow have been pe rf ormed by V e iierus , Schieber and c o worke rs 



Broerman et al 



(119991 ): 



Iddir et al. 



feood); 



Venerus et al. 



(|2Qo4); 



Venerus et al. 



Schieber et al 



1999); 



(120041 ): 



21 



Balasubramanian et al.l (120051 )] . and so me issues associated with generalizing the concept 



of heat capacity have been discussed by 



Hiitter et al 



fl2009h . 



Nonequihbrium thermodynamics guides the formulation of autonomous evolution equa- 
tions and offers the recipes for finding the thermodynamic building blocks on a less detailed 
level of description from a more detailed level of description. The thermodynamic structure is 
a geometric one and the passage between different levels of description hence is a m athemat- 



ical topic of structure preserving transformations [Ottinger and Struchtrupl (120071 )]. Ideally, 
the relationship between thermodynamics and mathematics should be even more intimate. 
For example, general mathematical results concerning a symptotic stability have already been 
obtained for metriplectic, or GENERIC, structures by iBirtea et al.l (120071 ). The criteria for 
thermodynamic admissibility, in particular, the existence of a nondecreasing entropy and 
a positive semidefinite friction matrix, should coincide with the criteria for mathematically 
proving the existence and uniqueness of solutions to the evolution equations. The existence 
of unique solutions is what we expect from meaningful equations. Thermodynamics shall be 
designed to guarantee that. 
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